Crystal nucleation of colloidal hard dumbbells 
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Using computer simulations we investigate the homogeneous crystal nucleation in suspensions 
of colloidal hard dumbbells. The free energy barriers are determined by Monte Carlo simulations 
using the umbrella sampling technique. We calculate the nucleation rates for the plastic crystal and 
the aperiodic crystal phase using the kinetic prefactor as determined from event driven molecular 
dynamics simulations. We find good agreement with the nucleation rates determined from sponta- 
neous nucleation events observed in event driven molecular dynamics simulations within error bars 
of one order of magnitude. We study the effect of aspect ratio of the dumbbells on the nucleation 
of plastic and aperiodic crystal phases and we also determine the structure of the critical nuclei. 
Moreover, we find that the nucleation of the aligned CP1 crystal phase is strongly suppressed by a 
high free energy barrier at low supersaturations and slow dynamics at high supersaturations. 

PACS numbers: 82.70.Dd, 64.60.Q-, 64.60.qe 



I. INTRODUCTION 

Recent breakthroughs in particle synthesis produced a 
spectacular variety of anisotropic building blocks^ Col- 
loidal particles with the shape of a dumbbell are one of 
the simplest anisotropic building blocks. Their unique 
morphologies lead to novel self-organized structures. For 
instance, it was found that magnetic colloidal dumbbells 
can form chain-like clusters with tunable chirality^ while 
novel crystal structures have been predicted recently for 
asymmetric dumbbell particles consisting of a tangent 
large and small hard sphere, which are atomic analogs of 
NaCl, CsCl, 7CuTi, CrB, and alrV, when we regard the 
two individual spheres of each dumbbell independently. 3 
Moreover, colloidal dumbbells also gain increasing scien- 
tific attention in recent years due to its potential use in 
photonic applications. It has been shown that dumbbells 
on a face-centered-cubic lattice where the spheres of the 
dumbbells form a diamond structure exhibit a complete 
band gap^£ while it is impossible to obtain a complete 
band gap in systems consisting of spherical particles. A 
very recent calculation showed that for midrange aspect 
ratios, both asymmetric and symmetric dumbbells have 
2-3 large band gaps in the inverted latticed Although 
these structures are not thermodynamically stable for 
hard dumbbells,— ~— it does show the promising poten- 
tial of anisotropic particles in photonic applications. 

New routes of synthesizing colloidal dumbbells make it 
easy to control the aspect ratioJ^ In addition, by adding 
salt to the solvent, the interactions between dumbbells 
can be tuned from long-ranged repulsive to hard interac- 
tions. Although hard dumbbells were originally modeled 
for simple non-spherical diatomic molecules, such as ni- 
trogen, they are also a natural model system for studying 
the self-assembly of colloidal dumbbells*^— The phase 
behavior of hard dumbbells has been extensively stud- 
ied by density functional theory^ and computer simu- 
lations^— The bulk phase diagram of hard dumbbells 
displays three types of stable crystal structures^— For 



small aspect ratio, the dumbbells form a plastic crystal 
phase at low densities. The freezing into a cubic plastic 
crystal phase in which the dumbbells are positioned on a 
face-centered-cubic lattice but are free to rotate, has been 
determined using Monte Carlo simulations^ These re- 
sults have been refined by Vega, Paras, and Monson, who 
showed that at higher densities the cubic plastic crystal 
phase transforms into an orientationally ordered crystal 
CP1 phase. Additionally, these authors showed that the 
fluid-cubic plastic crystal coexistence region terminates 
at L/a ~ 0.38, where L is the distance between the cen- 
ters of spheres and a is the diameter of the dumbbells. 
For longer dumbbells a fluid-CPl coexistence region was 
found, whereas the relative stability of the close-packed 
crystal structures CP1, CP2, and CP3, which only dif- 
fer in the way the hexagonally packed dumbbell layers 
are stacked remained undetermined as the free energies 
are very similar J£r— Moreover, these authors showed by 
making an estimate for the degeneracy contribution to 
the free energy that dumbbells with L/a = 1 may form 
an aperiodic crystal phased The stability of such an ape- 
riodic crystal structure in which both the orientations 
and positions of the particles are disordered, while the 
spheres of each dumbbell are located on the lattice posi- 
tions of a random-hexagonal-close-packed (rhcp) lattice, 
has been verified recently for L/a > 0.88^ In addition, 
it has been shown that the plastic crystal phase with the 
hexagonal-close-packed structure is more stable than the 
cubic plastic crystal for a large part of the stable plas- 
tic crystal region^ 3 Although, the bulk phase diagram is 
well-studied, the kinetic pathways of the fluid-solid phase 
transitions are still unknown, and only a few studies have 
been devoted to the crystal nucleation of anisotropic par- 
ticles^^ In the present work, we investigate the nucle- 
ation of the plastic crystal phase of hard dumbbells using 
computer simulations and study the effect of aspect ratio 
of the dumbbells on the resulting nucleation rates and 
the structure and size of the critical nuclei. Moreover, 
for longer dumbbells we investigate crystal nucleation of 
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the aperiodic crystal phase. First, we calculate the free 
energy barriers for nucleation using Monte Carlo (MC) 
simulations with the umbrella sampling technique, which 
are then combined with event driven molecular dynam- 
ics (EDMD) simulations to determine the kinetic prefac- 
tor and the nucleation rates. Additionally, we determine 
the nucleation rates from spontaneous nucleation events 
observed in EDMD simulations. We compare the nucle- 
ation rates and critical nuclei obtained from the umbrella 
sampling MC simulations with those from EDMD simu- 
lations. 

The remainder of this paper is organized as follows. In 
Sec. II, we describe the methodology including the model 
and simulation methods used. We present the results and 
discussions on the nucleation of three types of crystal 
phases in suspensions of hard dumbbells in Sec. III. We 
end with some discussions and conclude in Sec. IV. 



II. METHODOLOGY 

We consider a system of hard dumbbells consisting of 
two overlapping hard spheres with diameter a with the 
centers separated by a distance L. We define the aspect 
ratio as L* = L/a, such that the model reduces to hard 
spheres for L* = and to tangent spheres for L* = 1. 
We study crystal nucleation of hard dumbbells for < 
L* < 1. We focus on the nucleation of the plastic crystal 
phase (0 < L* < 0.4) and the aperiodic crystal phase 

(o.88 < l* < i)*ia-ia 



A. Order parameter 

In order to study the nucleation of the crystal phase, 
we require a cluster criterion that identifies the crystalline 
clusters in a metastable fluid. In this work, we employ the 
order parameter based on the local bond order parameter 
analysis of Steinhardt et al.S^ We define for every par- 
ticle i, a 21 + 1-dimensional complex vector q;(i) given 
by 
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where Nb(i) is the total number of neighboring particles 
of particle i, and T; m (fy ) is the spherical harmonics for 
the normalized direction vector f y between particle i and 
j, I is a free integer parameter, and m is an integer that 
runs from m — — I to m = +1. Neighbors of particle i are 
defined as those particles which lie within a given cutoff 
radius r c from particle i. In order to determine the cor- 
relation between the local environments of particle i and 
j, we define the rotationally invariant function di(i,j) 



dl(i,j) = qi m (i) ■ q~L(j) 



(2) 



where qi m (i) = gim(i)/yJ2 m =-l \lim{i)\ 2 and the aster- 
isk is the complex conjugate ? 2 If di(i,j) > d c , the bond 
between particle (sphere) i and j is regarded to be solid- 
like or connected, where d c is the dot-product cutoff. We 
identify a particle (sphere) as solid-like when it has at 
least £ c solid-like bonds. We have chosen the symmetry 
index ( = 6 as the particles (spheres) display hexago- 
nal order in the plastic crystal and the aperiodic crystal 
phase. We have chosen r c = 1.3a, d c = 0.7, and £ c = 6 
in our simulations. It has been shown recently that the 
choice of order parameter (r c , d c , and £ c ) does not af- 
fect the resulting nucleation rate if it is not too restric- 



tive 



20,23 



To analyze the structure of the critical nuclei, we use 
the averaged local bond order parameter q l and Wi pro- 
posed by Lechner and DellagOfS 4 which allows us to iden- 
tify each particle as fcc-like or hep-like, provided the 
number of neighboring particles Nb(i) > 10: 
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21 + 



(3) 
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wi{i) = 



where 



mi+m2+m 3 =0 



3/2 



(4) 
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m= — l 



The sum from k = to Nb(i) runs over all neighbors 
of particle (sphere) i plus the particle (sphere) i itself. 
While qi m (i) takes into account the structure of the first 
shell around particle i, the averaged q lm {i), contains also 
the information of the structure of the second shell, which 
increases the accuracy of the crystal structure determi- 
nation. In order to distinguish fcc-like and hep-like parti- 
cles, we employ q 4 and w\, as the order parameter distri- 
butions of pure fee and hep phases of Lennard- Jones and 
Gaussian core systems are well separated in the q 4 — W4 
planed 



B. Umbrella sampling 

The Gibbs free energy AG(n) for the formation of a 
crystalline cluster of size n is given by AG(n)/fcsT = 
const — ln[P(n)], where P(n) is the probability distribu- 
tion function of finding a cluster of size n, fcs is Boltz- 
mann's constant, and T the temperature. As nucleation 
is a rare event and the probability to find a spontaneous 
nucleation event is very small in a brute force simulation 
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within a reasonable time, one has to resort to special- 
ized simulation techniques such as forward flux sampling, 
umbrella sampling or transition path sampling. Here, we 
employ the method developed by Frenkel and cowork- 
ers^ to calculate the free energy of the largest cluster. 
In this method, the sampling is biased towards configu- 
rations that contain clusters with a certain size. To this 
end, we introduce a biasing potential oj(r N ), which is a 
harmonic function of the cluster size n: 



PLu(r N ) = h[n(r N )-n ] 2 



(6) 



where n(r N ) is the size of largest cluster and no is the 
center of the umbrella sampling window whose width de- 
pends on k. In this work we set k = 0.2. By increasing 
the value of no, we increase the size of the largest crys- 
talline cluster in our system, which enables us to cross 
the nucleation barrier. If we define the average number of 
crystalline clusters with n particles by (N n ), one can cal- 
culate the probability distribution P(n) — {N n )/N from 
which we can determine the Gibbs free energy AG(n). 
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FIG. 1: Gibbs free energy AG(n) /ksT as a function of cluster 
size n for the nucleation of the plastic crystal phase of hard 
dumbbells with various aspect ratios L* — L/a as displayed 
and supersaturation /3|A/x| = 0.34 (filled symbols) and 0.54 
(open symbols). 



C. Event driven molecular dynamics simulations 

Since the potential between particles in systems of 
hard dumbbells is discontinuous, the pair interactions 
only change when particles collide. The particles per- 
form elastic collisions when they encounter each other. 
We numerically identify and handle these collisions by 
using an EDMD simulation! 26 ! 27 

Using MD simulations to determine the nucleation rate 
is straightforward. Starting with an equilibrated fluid 
configuration, an MD simulation is used to evolve the 
system until the largest cluster in the system exceeds the 
critical nucleus size. Then the nucleation rate is given by 

where (t) is the averaged waiting time of forming a crit- 
ical nucleus in a system of volume V. 

III. RESULTS AND DISCUSSIONS 

In this section, we present the results on the nucleation 
of the plastic crystal, the aperiodic crystal and the CP1 
crystal phase in suspensions of hard dumbbells. 

A. Nucleation of the plastic crystal phase 

We first investigate the nucleation of the plastic crys- 
tal phase of hard dumbbells. Monte Carlo simulations 
with the umbrella sampling technique are performed on 
hard-dumbbell fluids with L* — 0,0.15 and 0.3 at su- 
persaturation /3|A/x| = 0.34 and with L* = 0,0.15 and 
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FIG. 2: Mean square displacement (Ar 2 (t)} as a function of 
time t/r in a fluid of hard dumbbells with L* = 0.3 at P* — 30 
(for p\Afi\ = 0.47). 



0.2 for (3\Au\ = 0.54 with /3 = l/k B T. We have chosen 
a shorter aspect ratio for the highest supersaturation as 
the plastic crystal phase for dumbbells with L* = 0.3 be- 
comes metastable with respect to the aligned CP1 phase 
for P* = Pa 3 /k B T > 30, i.e., p\Afi\ > 0.47. The Gibbs 
free energy (3AG{n) as a function of cluster size n is 
shown in Fig. Q] We clearly observe that at low super- 
saturation, i.e. /3|A/x| = 0.34, the heights of the free 
energy barriers increase slightly (~ 8%) with aspect ra- 
tio. More specifically, (3 AG* =42.9 ± 0.3, 44.5 ±1.1, and 
45.2±2 for L* = 0,0.15 and 0.3, respectively. Accord- 
ing to classical nucleation theory (CNT), the nucleation 
barrier for a spherical nucleus with radius R is given by 
AG(R) = 4n-fR 2 - 4tt| A[i\p s R 3 /3 with 7 the interfa- 
cial tension, \Afi\ the chemical potential difference be- 
tween the solid and fluid phase, and p s the bulk den- 
sity of the solid phase. CNT predicts a nucleation bar- 
rier height AG* = (167T/3)7 3 /( j o s |A/j,|) 2 and a critical 
radius R* — 2j / p s \Afj,\. The small increase in barrier 
height with aspect ratio can be explained by the small 
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increase in the crystal-melt interfacial tensions that have 
been determined recently for the crystal planes (100), 
(110), (111) using nonequilibrium work measurements 
with a cleaving procedure in MC simulations^ For a 
spherical cluster, the surface tension is expected to be an 
average over the crystal planes, i.e., (3jd 2 — 0.58,0.57, 
and 0.60, for L* = 0,0.15 and 0.3, respectively, where 
d 3 = cr 3 (l + 3/2L* - 1/5L* 3 ). Another work by David- 
chack et al. found a slightly lower value for the averaged 
interfacial tension of hard spheres, i.e, f3^d 2 = 0.559. 29 
Using /?AG* = 42.9 and the more precise value for the 
surface tension f3jd 2 — 0.559, and the values for f3jd 2 
and the bulk density p s for varying L* presented in Ta- 
ble U CNT predicts a slightly lar ger increase in barrier 
height upon increasing L* , i.e., (3 AG* = 45.6 and 50.49, 
for L* = 0.15 and 0.3, respectively. However, when the 
supersaturation is increased to /3|A/i| = 0.54, we find 
a decrease in barrier height upon increasing the aspect 
ratio as shown in Fig. [1] and Table HI which cannot be ex- 
plained by CNT. Apparently, the pressure dependence of 
the surface tension is different for dumbbells with various 
aspect ratios. 

The nucleation barriers obtained from umbrella sam- 
pling MC simulations can also be used to determine the 
nucleation rates as given by— 

7 = /texp(-/3AG*) (8) 

where k is the kinetic prcfactor given by n — 
Plfn* \/\AG" (n*)\/2TrkBT , p\ is the number density of 
particles in the fluid phase, /„* the rate at which par- 
ticles are attached to the critical nucleus, AG (n*) is 
the second derivative on the top of the Gibbs free energy 
barrier. The attachment rate can be calculated from the 
mean square deviation of the cluster size at the top of 
the free energy barrier by 

_l ([n(t)-n(0)} 2 ) 
h " ~ 2 1 (9) 

where n(t) is the cluster size at time t. The mean square 
deviation of the cluster size can be determined from 
EDMD simulations starting from configurations at the 
top of the free energy barriers. Using the results for the 
attachment rates and the nucleation barriers obtained 
from umbrella sampling MC simulations, we can deter- 
mine the nucleation rates, which we compare with those 
obtained directly from spontaneous nucleation events in 
EDMD simulations. We observed a large variance in the 
attachment rates calculated for different nuclei. We used 
10 independent configurations on the top of the barrier 
and followed 10 trajectories for each of them to deter- 
mine the attachment rates. Taking into account the sta- 
tistical errors in the free energy barriers and attachment 
rates, we estimate that the error in the resulting nucle- 
ation rates is one order of magnitude. In order to exclude 
the effect of dynamics, we compare the nucleation rates 
for the plastic crystal phase in long-time diffusion times, 
i.e. tl = a 2 /6Di with Di the long-time diffusion coef- 
ficient. We calculate D\ by measuring the mean square 



displacement at supersaturation /3|A/z| = 0.34 and 0.54 
as shown in Table U for various aspect ratios. We clearly 
observe that the dynamics becomes slower for increasing 
aspect ratio L*, resulting in long-time diffusion coeffi- 
cients Dir/a 2 = 0.012,0.01,0.0023 for L* = 0,0.15 and 
0.3 at /3|A/i| = 0.34 with r = a^fm/k B T. At higher su- 
persaturation /3|A/i| = 0.54, we find even smaller values 
for i-e., Dir/a 2 = 0.0078, 0.006, 0.003 for L* = 0,0.15 
and 0.2, respectively. 

The resulting nucleation rates in units of the long-time 
diffusion coefficient are shown in Table |U We wish to 
make a few remarks here. First, the nucleation rates 
obtained from spontaneous nucleation events observed 
in EDMD simulations agree well with the ones obtained 
from umbrella sampling MC simulations within error bars 
of one order of magnitude, which means that the nucle- 
ation results obtained from the umbrella sampling MC 
simulations are reliable. Secondly, we clearly observe 
that the nucleation rates for the different aspect ratios 
ranging from L* — to 0.3 are remarkably similar as the 
differences are within the errorbars for both supersatu- 
rations. 

Finally, we made an attempt to study spontaneous nu- 
cleation of dumbbells with L* = 0.3 at supersaturation 
/3|A/i| — 0.54 using event-driven MD simulations. As 
already mentioned above, the plastic crystal phase for 
dumbbells with L* = 0.3 becomes metastable with re- 
spect to the aligned CP1 phase for P* = Pa 3 /k B T > 
30, i.e., /3|A/i| > 0.47. Hence, we would expect to 
find the nucleation of the CP1 phase here. How- 
ever, we find that the nucleation is severely hampered 
due to slow dynamics, which can be appreciated from 
Fig. [21 where we plot the mean square displacement for 
/3|A/i| = 0.47. The resulting long-time diffusion coeffi- 
cient Di = 1.72 x 10~ 4 ct 2 /t is at least one order of mag- 
nitude smaller than the long-time diffusion coefficients at 
/3| A^| = 0.54, where we observed spontaneous nucleation 
for L* = 0,0.5, and 0.2. 

In umbrella sampling MC simulations, we can "fix" the 
simulations at the top of the nucleation barrier which al- 
lows us to study the properties of the critical nuclei. We 
investigate the effect of the particle anisotropy on the 
structure of the critical nuclei using the order param- 
eters g 4 and u>4 as defined above. At supersaturation 
/3|A/i| = 0.34, the size of the critical nuclei is n ~ 250 
which is sufficiently large to determine the crystal struc- 
ture of the nuclei. For each dumbbell, we calculate the 
averaged local bond order parameter q 4 and W4, pro- 
vided the particle has iVb(z) > 10 neighbors. The dis- 
tribution of particles in the critical nuclei are presented 
as scatter plots in the g 4 — 1ZJ4 plane along with those 
for pure fee and hep plastic crystal phases of dumbbells 
with L* = 0,0.15, and 0.3, at corresponding pressures. 
From Fig. [31 we clearly observe that the critical nu- 
clei for L* = and 0.15, contains predominantly fcc-like 
rather than hep-like particles. In order to distinguish the 
fcc-like and hep-like particles more quantitatively, we di- 
vide the g 4 — W4 plane by a straight line in such a way 



5 



IS o 







1 

■ FCC 

■ HCP 

critical nucleus 







IS 



0.05 O.I _ 0.15 0.2 



(b) 




1.9q 4 -w 4 = 0.35/ 



FCC 
HCP 

critical nucleus 



IS 



0.05 O.I _ 0.15 0.2 

% 




FIG. 3: Distribution of particles in the critical nuclei for the plastic crystal nucleation in hard dumbbell systems with L* = 
(a), 0.15 (b), and 0.3 (c) as obtained from umbrella sampling MC simulations at a supersaturation /3|A/i| = 0.34 in the q 4 — UJ4 
plane compared with those for pure fee and hep plastic crystal phases with corresponding pressures. The dashed lines are used 
to distinguish the fee-like and hep-like particles, and the formulas are next to them. 






FIG. 4: Typical configurations of critical nuclei for the plastic crystal nucleation of hard dumbbells with aspect ratios L* = 
(a), 0.15 (b), and 0.3 (c) at supersaturation |A/i| = 0.34fc_gT. The red (dark grey) particles are fee-like, the blue particles are 
hep-like particles, while the light blue (light grey) particles are undetermined. 



that the particle distributions for the pure fee and hep 
plastic crystal phases are maximally separated. We plot 
the criteria to distinguish fcc-like and hep-like particles 
as dashed straight lines in Fig. [3] with the correspond- 
ing formula. We note, however, that the criteria seem 
to be arbitrarily chosen, but the identification of fcc-like 
and hep-like particles for typical nuclei seems to be less 
sensitive on the precise details of these criteria. Typical 
snapshots of the critical nuclei for L* = 0,0.15, and 0.3 
are shown in Fig. 21 where the color-coding denotes the 
identity (fcc-like, hep-like or undetermined) of the parti- 
cle using these criteria. As we did not calculate the aver- 
aged local bond order parameter g 4 and W4 for particles 
with Nb(i) < 10 neighbors, the identity of these particles 
remains undetermined. We clearly observe that the crit- 
ical nuclei for L* = 0,0.15 contains mainly fcc-like parti- 
cles. The particle distributions becomes broader for the 
pure fee and hep plastic crystal phases upon increasing L* 
and consequently it becomes more difficult to distinguish 
fcc-like and hep-like particles. However, the fraction of 
hep-like particles seems to increase with increasing par- 
ticle elongation. This agrees with the results from free 
energy calculations of hard dumbbell systems, where it 
has been shown that the hep plastic crystal phase is more 
stable than the one with an fee structure at L* > 0.15. 13 
It is worth noting here that recent nucleation studies of 
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FIG. 5: Gibbs free energy AG(n) as a function of number of 
spheres n in the largest cluster for the nucleation of the (ape- 
riodic) crystal phase of hard dumbbells and of hard spheres. 



hard spheres showed that the critical nuclei contain ap- 
proximately 80% fcc-like particles M° As the free energy 
difference per particle between bulk fee and hep phases is 
only about 0.001 ksT at melting, the predominance for 
fcc-like particles is attributed to surface effects. 
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TABLE I: Nucleation rates la /6D; for the nucleation of the plastic crystal phase in systems of hard dumbbells with elongation 
L* , at pressure Pa 3 /IcbT ', and supersaturation /3|A/ij. p s d 3 is the number density of dumbbells in the solid phase, /3AG(n*) 
is the barrier height, and \/3AG (n*)\ is the second derivative of the Gibbs free energy at the critical nucleus size n* , i.e., the 
number of dumbbells in the critical cluster. / n * /6Di is the attachment rate in units of the long-time diffusion coefficient Di. 
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FIG. 6: Equation of state (EOS, i.e., /3Pa vs the number 
density of spheres p sph a 3 for a system of hard spheres and 
hard dumbbells with L* — 1.0. For the fluid and solid phase 
of hard spheres, the Carnahan-Starling 30 and Speedy 31 EOS 
are plotted. The EOS of hard dumbbells for the fluid phase is 
obtained from Ref. 1321 . The dashed vertical line denotes the 
bulk coexistence pressure of hard dumbbells with L* — 1.0. 



B. Nucleation of the aperiodic crystal phase 

For more elongated dumbbells, i.e. L* > 0.88, the ori- 
cntationally disordered aperiodic crystal phase becomes 
stable^""— in which the individual spheres of the dumb- 
bells are on a random hep lattice whereas the orienta- 
tions of the dumbbells are random. In this section, we 
investigate the nucleation of the aperiodic crystal phase 
of hard dumbbells with different aspect ratios. We per- 
form Monte Carlo simulations using the umbrella sam- 
pling technique to determine the Gibbs free energy as a 
function of cluster size for hard dumbbells with L* = 1.0 
and supersaturation P* = 16 and 17. The order pa- 
rameter that is employed here in the umbrella sampling 
technique is equal to the number of spheres n (and thus 
not the number of dumbbells) in the largest crystalline 
cluster in the system. Thus, we check for each individual 
sphere whether or not it belongs to the largest crystalline 



cluster, and as a consequence, the whole dumbbell can be 
part of the largest cluster or only one sphere of the dumb- 
bell can belong to the cluster, or the whole dumbbell is 
regarded to be fluid-like. Consequently, it is convenient 
to introduce a bulk chemical potential per sphere, which 
equals 0.5 times the bulk chemical potential per dumb- 
bell pi sph — fi/2. We compare the results with those for 
hard spheres at the same pressure in Fig. [5] Since the 
bulk pressure for the solid-fluid transition of hard dumb- 
bells with L* — 1 is remarkably close to that of hard 
spheres f3P coex a 3 = 11.8^""— one might naively expect 
that the nucleation barriers should be compared at the 
same dimensionless pressure. However, we observe that 
at the same pressure, the nucleation barrier for the ape- 
riodic crystal phase of hard dumbbells is slightly higher 
than that of hard spheres. CNT predicts that the barrier 
height is given by AG* = (167r/3)7 3 /(p^|A^ s P /l |) 2 , and 
hence a difference in barrier height should be due to a dif- 
ference in the interfacial tension 7, the density of spheres 
in the solid phase p s s ph , or in |A// sp ' l |. As the reduced 
density of spheres p sph a 3 in the aperiodic crystal phase 
is very close to that of a solid phase of hard spheres at 
P* = 16 and 17, and the interfacial tensions f3"fcr 2 are 
also expected to be very similar, the difference in barrier 
height can only be caused by a difference in |A/i sp ' l |. We 
therefore calculated more accurately the bulk chemical 
potential difference per sphere between the solid and the 
fluid phase using 



|A/i 



sph I 



1 



1 



sph 



sph 



dp 



(10) 



where pf and p s s ph are the density of spheres in the 
liquid and solid phase. In Fig. [6J we plot the equation of 
state for the fluid and solid phase of hard spheres from 
Ref. [3fj||3ll along with the equation of state for the fluid 
phase of hard dumbbells for L* = 1 from Ref. [HI- In 
addition, we determined the equation of state for the solid 
phase using EDMD simulations. Using these results and 
Eq. rrjJl we indeed find that the supersaturation f3\Ap sph \ 
per sphere is ~ 2.3% smaller for hard dumbbells than for 
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FIG. 7: Number of spheres in the largest cluster n max as a 
function of time t/r for a typical trajectory obtained from 
EDMD simulations for the aperiodic crystal nucleation of 
hard dumbbells with L* = 1.0, N = 16000 and P* = 17. 



FIG. 9: Mean square deviation of the cluster size < [n(t) — 
n(0)] 2 > as a function of time t/r for hard dumbbells with 
L* = 0.95,0.97, and 1.0 at supersaturation /3|A^ sph | = 0.43. 
The resulting attachment rates /„* are listed in units of r _1 . 
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FIG. 8: Gibbs free energy AG(n) as a function of the number 
of spheres n in the largest crystalline cluster for the aperiodic 
crystal nucleation of hard dumbbells with L* = 0.95, 0.97, 
and 1 at supersaturation f3\A(i sph \ = 0.43. 



hard spheres, resulting in an increase in barrier height of 
~ 5%, which perfectly matches our results. We conclude 
that the difference in the height of the nucleation barrier 
between the aperiodic crystal phase of dumbbells with 
L* = 1.0 and the hard-sphere crystal is mostly due to 
the difference in \Ap, sph \. 

Moreover, we also performed EDMD simulations for 
the spontaneous nucleation of the aperiodic crystal phase 
of dumbbells at P* = 17 in system of N = 16000 hard 
dumbbells. The number of spheres in the biggest cluster 
as a function of time from a typical MD simulation is 
shown in Fig. [7J We find that the size of critical nuclei 
in spontaneous nucleation is around 100 spheres which 
agrees well with the result obtained from umbrella sam- 
pling MC simulations shown in Fig. [5] The nucleation 
rate obtained from spontaneous nucleation events ob- 
served in MD simulations is Ia 5 /6Di = 7.3 x 10 _8±1 
which agrees very well with the rate obtained from um- 
brella sampling MC simulations, Ia 5 /6Di = 2.8xl0" 8±1 , 




FIG. 10: Nucleation rate Ia 5 /6Di for the aperiodic crystal 
phase as a function of the aspect ratio L* of hard dumbbells 
at supersaturation /3\A{i sph \ = 0.43. 



within the error bars of one order of magnitude. 

Furthermore, we study the effect of aspect ratio on 
the nucleation of the aperiodic crystal phase, and the 
free energy barriers for hard dumbbells with aspect ratios 
L* = 0.95,0.97, and 1.0 at supersaturation (3\A^ s P h \ = 

0. 43. We plot AG(n) as a function of cluster size n, 

1. e., the number of spheres in the cluster, in Fig. [8] 
We observe that at the same supersaturation the bar- 
rier height decreases upon decreasing the elongation of 
the dumbbells. According to classical nucleation theory, 
AG* oc ^/{pf h \A^P h \) 2 : where A^ h is the supersat- 
uration per sphere with p s s ph the bulk density of spheres 
in the solid phase. As shown in Table [TTl p s s ph a 3 is very 
similar for L* — 0.95,0.97, and 1.0, and we argue that 
the interfacial tension of the aperiodic crystal decreases 
upon decreasing the elongation of the dumbbells. In or- 
der to calculate the nucleation rates, we perform EDMD 
simulations starting from configurations on the top of the 
free energy barriers. We plot the mean square deviation 
of the cluster size as a function of time in Fig. [SJ We find 
that the attachment rate decreases significantly as the 
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FIG. 11: Distribution of the spheres in the critical nuclei as 
obtained from umbrella sampling MC simulations in q A — W4 
plane in systems of hard dumbbells with V — 1.0 at su- 
persaturation /3\Ap sph \ = 0.43. Inset: Typical configuration 
of a critical nucleus. Red denotes fcc-like spheres and blue 
denotes hep-like spheres while the light blue are the undeter- 
mined ones. 



anisotropy of the dumbbells decreases. The resulting nu- 
cleation rates in units of the long time diffusion coefficient 
are shown in Fig. [TO] We clearly observe that at fixed su- 
persaturation the nucleation rate increases with decreas- 
ing dumbbell elongation. However, in the phase diagram 
of hard dumbbells^ - — the pressure range where the ape- 
riodic crystal phase is thermodynamically stable shrinks 
significantly when the aspect ratio decreases. As a result, 
it is not possible to increase the supersaturation further 
for shorter dumbbells, although the nucleation rates are 
already much higher for shorter ones than for longer ones 
at the same supersaturation. 

Additionally, we also study the structure of the critical 
nuclei by calculating the averaged local bond order pa- 
rameter q 4 and W4, provided the sphere has Nb(i) > 10 
neighbors. The distribution of spheres in the critical nu- 
clei are presented as scatter plots in the q 4 — W4 plane 
in Fig. Qj] for L* = 1.0. We observe only a few spheres 
with W4 > and q 4 < 0.1, as most of the spheres are 
in the area of W4 < and q 4 > 0.1, which is very simi- 
lar to the scatter plots for hard spheres shown in Fig. |3^. 
Consequently, the critical nucleus of the aperiodic crystal 
phase of hard dumbbells contains also more fcc-like than 
hep-like particles, similar to the critical nuclei observed 
in hard-sphere nucleation^ A typical configuration of a 
critical nucleus is shown in the inset of Fig. [TTJ where the 
spheres are considered to be fcc-like if W4 < 0. 



C. Slow dynamics of hard dumbbells 

The phase diagram of hard dumbbells shows a sta- 
ble aligned CP1 crystal phase at infinite pressure for all 
aspect ratios of the dumbbells, and a fluid-CPl coexis- 
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FIG. 12: Estimated height of the Gibbs free energy barrier 
AG* jksT obtained from classical nucleation theory (solid 
line) and the packing fraction r\ in the supersaturated fluid 
phase' 12 (dashed line) as a function of pressure P* for the nu- 
cleation of the CP1 phase of hard dumbbells with V = 0.4. 



tence region for 0.4 < L* < 0.8pi2r— The surface ten- 
sion for the fluid-CPl interface of hard dumbbells with 
L* = 0.4 is f]ja 2 ~ 1.8i2£ The height of free energy bar- 
rier is given by AG* = 16nj 3 /3{p s \Afi\) 2 in CNT. If we 
assume that the interfacial tension does not change sig- 
nificantly with increasing pressure, we can estimate the 
free energy barrier height as a function of pressure by 
integrating the Gibbs-Duhem equation to obtain |A/x|. 
The barrier height AG* and the packing fraction 77 for 
the fluid phase are shown in Fig. [T^] as a function of the 
pressure P* . We find that the barrier height AG* is 
extremely high, and only becomes less than 50/cbT for 
P* > 45, corresponding to a packing fraction of the fluid 
phase rj > 0.67. However, if the interfacial tension in- 
creases with increasing pressure as shown in Ref. l33l . the 
"actual" height of free energy barrier can become even 
higher. As a consequence nucleation of the CP1 crystal 
phase is an extremely rare event. 

Additionally, we calculate mean square displacements 
(Ar 2 (t)) and the second-order orientational correlator 
L2(t) — (■f , 2[cos(0(t))]) for a metastable fluid of hard 
dumbbells with L* = 0.4, 0.5, and 0.8 at supersatura- 
tion /3|A/i| = 1.0 as shown in Fig. [131 We find that at a 
supersaturation /?|A/i| = 1.0, where the barrier height is 
still very high, AG*/k B T - 170 for L* = 0.4, the long- 
time diffusion coefficients D\ ~ 10~ 4 <7 2 /V obtained from 
(Ar 2 (i)) is extremely small, see Table III), whereas L 2 (t) 
exhibits slow relaxation. Our findings are consistent with 
predictions obtained from mode-coupling theory for a 
liquid-glass transition, in which the structural arrest is 
due to steric hindrance for both translational and reori- 
cntational motioni^~— Moreover, mode-coupling theory 
predicts that the steric hindrance for reorientations be- 
comes stronger with increasing elongation, which is con- 
sistent with our results for L 2 (t) in Fig. [T5t^~— . In- 
creasing the supersaturation will lower AG* , but Di will 
decrease as well, while at lower supersaturation the bar- 
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TABLE II: Nucleation rates la /6-D; for the nucleation of the aperiodic crystal phase in systems of hard dumbbells with 
elongation L* , at pressure Pa 3 /knT, and supersaturation per sphere P\Ap, svh \. p B s ph a 3 is the number density of spheres in the 
solid phase, /3AG(n*) is the barrier height, and \/3AG (n*)\ is the second derivative of the Gibbs free energy at the critical 
nucleus size n* , i.e., the number of spheres in the critical cluster. / n * /6D; is the attachment rate in units of the long-time 
diffusion coefficient Di. 



L* 


Pa 3 /k B T 


f3\Afi aph \ 


P7 h <r 3 


n* 


/3AG(n*) 


\j3AG"(n*)\ 


/n-/6A 


Ar/a 2 


7a 5 /6A (US) 


Ia 5 /6Di (MD) 


1 


17 


0.53 


1.170 


115 


21.4 ±0.4 


1.2 x 10~ 3 


2813 


0.0026 


2.0 x 10~ 8±1 


7.3 x 10" 8±1 


1 


16 


0.43 


1.158 


170 


29.5 ±0.6 


9.4 x 10 -4 


5556 


0.0036 


1.1 x 10 _11±1 




0.97 


18 


0.43 


1.171 


140 


25.3 ±0.9 


8.4 x 10~ 4 


5228 


0.0022 


6.6 x 10~ 10±1 




0.95 


20 


0.43 


1.182 


100 


19.9 ±0.7 


3.0 x 10 -3 


2273 


0.0011 


1.2 x 10 _7±1 





TABLE III: Long -time d iffusion coefficients Di in units of 
cr 2 /r with r = a^/m/kBT for hard dumbbells with elongation 
L* at pressure P* , packing fraction -q, and supersaturation 
/3|A M | = 1.0. 



L* 


P* 


■n 


Dit/g 2 


0.4 


34.5 


0.64 


1.02 x 10~ 4 


0.5 


31.2 


0.63 


2.47 x 10 -4 


0.8 


24.8 


0.61 


2.78 x 10 -4 



rier height will only increase. As a result, the nucleation 
of CP1 phase of hard dumbbells is severely hindered by 
a high free energy barrier at low supersaturations and 
slow dynamics at high supersaturations, which explains 
why the CP1 phase of colloidal hard dumbbells has never 
been observed in experiments^ or in direct simulations. 
It is worth noting that the phase diagram might also 
display (meta)stable CP2 and CP3 close-packed crystal 
structures^ which only differ in the way the hexago- 
nally packed dumbbell layers are stacked. As the free 
energy difference for the three close-packed structures 
is extremely small, we expect the surface tensions and 
the nucleation barrier height to be very similar. Hence, 
we expect that also the nucleation of the CP2 and CP3 
phases are hindered by either a high free energy barrier 
or slow dynamics. 



IV. CONCLUSIONS 

In conclusion, we investigated the homogeneous nu- 
cleation of the plastic crystal, aperiodic crystal and 
CP1 crystal phase of hard dumbbells using computer 
simulations. Hard dumbbells serve as a model sys- 
tem for colloidal dumbbells for which the self-assembly 
is mainly determined by excluded volume interactions. 
For charged colloidal dumbbells or diatomic molecules, 
screened Coulombic interactions and Van der Waals in- 
teractions may significantly change the kinetic pathways 
for nucleation. For instance, crystal nucleation of hard 
rods proceeds via multi-layered crystalline nuclei whereas 
attractive depletion interactions between the rods in a 
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FIG. 13: Mean square displacements (Ar 2 (t)) and L2(t) = 
(P2[cos(8(t))]) as a function of time for a fluid of hard dumb- 
bells with an aspect ratio L* = 0.4, 0.5, and 0.8 at supersat- 
uration /3jA/^| = 1.0. 

polymer solutions favor the nucleation of single-layered 
nuclei. 20 ' 39 For the nucleation of the plastic crystal phase 
of hard dumbbells, we found that at low supersatura- 
tions the free energy barriers increases slightly with in- 
creasing dumbbell anisotropy, which can be explained by 
a small increase in surface tension for more anisotropic 
dumbbells. 28 When the supersaturation increases, the 
barrier height decreases with increasing dumbbell aspect 
ratio, which can only be explained by a different pressure- 
dependence of the interfacial tension for hard dumbbells 
with different aspect ratios. Although the nucleation 
rate for the plastic crystal phase does not vary much 
with aspect ratio, the dynamics do decrease significantly. 
We also carried out EDMD simulations and compared 
the nucleation rates obtained from spontaneous nucle- 
ation events with those obtained from the umbrella sam- 
pling Monte Carlo simulations, and found good agree- 
ment within the error bars of one order of magnitude. 
Additionally, we investigated the structure of the criti- 
cal nuclei of the plastic crystal phase of hard dumbbells 
with various aspect ratios. We found that the nuclei of 
the plastic crystal tend to include more fcc-like particles 
rather than hep- like ones, which is similar to the critical 
nuclei of hard spheres^ However, the amount of hep-like 
particles increases with increasing dumbbell aspect ratio, 
which agrees with the free energy calculations^ where it 
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has been shown that the hep structure is more stable 
than fee structure for L* > 0.15. 

Moreover, we also studied the nucleation of the ape- 
riodic crystal phase of hard dumbbells, and our results 
showed that at the same pressure, the nucleation barrier 
of the aperiodic crystal phase of hard dumbbells with 
L* = 1.0 is slightly higher than that of hard spheres 
which is mostly due to a small difference in supersatu- 
ration (3\A[i sph \. We also performed EDMD simulations 
for the spontaneous nucleation of the aperiodic crystal 
from hard-dumbbell fluid phase, and we found that the 
nucleation rate obtained from spontaneous nucleation 
agrees very well with the one obtained from umbrella 
sampling MC simulations. Furthermore, we studied the 
effect of aspect ratio on the nucleation of the aperiodic 
crystal phase, and found that at the same supersatura- 
tion, the nucleation rate in units of long-time diffusion 
coefficients increases for shorter hard dumbbells. How- 
ever, when the aspect ratio of dumbbells decreases, the 
pressure range where the aperiodic crystal phase is sta- 
ble becomes smaller. Additionally, we also found that 
the structure of the critical nuclei of the aperiodic crys- 



tal phase formed by hard dumbbells with L* — 1.0 is 
very similar to that of hard spheres which tend to have 
more fcc-like particles rather than hep-like ones. 

We estimated the height of the free energy barrier for 
the nucleation of the CP1 crystal phase of hard dumbbells 
according to classical nucleation theory, which turns out 
to be extremely high in the normal pressure range due 
to a high interfacial tension. Furthermore, we calculated 
the long-time diffusion coefficients for hard dumbbells at 
a moderate supersaturation, i.e. /3|A/i| = 1.0, which ap- 
pears to be very small. As a result, we conclude that 
the high free energy barrier as well as the slow dynamics 
suppress significantly the nucleation of CP1 phase. 
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